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ABSTRACT: Anapole moments are parity-odd, time-reversal-even moments of the El projec- 
tion of the electromagnetic current. Although it was recognized, soon after the discovery of parity 
violation in the weak interaction, that elementary particles and composite systems like nuclei 
must have anapole moments, it proved difficult to isolate this weak radiative correction. The first 
successful measurement, an extraction of the nuclear anapole moment of 133 Cs from the hyperfine 
dependence of the atomic parity violation, was obtained only recently. An important anapole 
^ | moment bound in Tl also exists. We discuss these measurements and their significance as tests 

of the hadronic weak interaction, focusing on the mechanisms that operate within the nucleus to 
generate the anapole moment. The atomic results place new constraints on weak meson- nucleon 
couplings, ones we compare to existing bounds from a variety of p — p and nuclear tests of parity 
nonconservation . 
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1 INTRODUCTION 



Until 1957 physicists assumed that the fundamental laws of nature did not distin- 
guish between left and right. However in that year, following a suggestion by Lee 
and Yang (||), experimenters discovered that the weak force governing processes 
such as muon decay and nuclear beta decay violated mirror symmetry maximally 
(to the accuracy of the measurements) (|2|). Very soon afterwards Vaks and Zel- 
dovich @ noted that weak interactions would then modify the electromagnetic 
couplings to elementary particles (as well as composite systems), allowing them 
to have a new, parity-violating moments, called anapole moments, in addition to 
their familiar parity-conserving ones (e.g., the charge and magnetic moments). 

This new moment has a number of curious properties. It vanishes when probed 
by real photons, and thus must be tested in processes where a virtual photon is 
exchanged. Thus, for example, the anapole moment of a nucleus can be probed in 
electron scattering and can influence the energies of bound electrons in an atom, 
but cannot be measured through direct interactions with an electric field (unlike 
magnetic moment interactions in a static B field). The resulting electron- nucleus 
interaction in an atom is pointlike, thus mimicking the short-range tree-level weak 
interaction induced by Z° exchange between atomic electrons and the nucleus. The 
atomic cloud feels the nuclear anapole moment only to the extent that the orbiting 
electrons penetrate the nucleus. The anapole moment is an electric dipole coupling 
that is nuclear spin-dependent: it is this spin dependence, as we will see, that 
allows anapole effects to be separated from tree-level weak interactions. Finally, 
the anapole moment is in general one of a larger class of weak radiative corrections. 
The anapole interaction - one diagram is given in Fig. la - is thus accompanied 
by other radiative corrections that do not correspond to virtual photon exchange, 
such as Fig. lb. It is the sum of all such diagrams that contribute to physical 
observables. It follows that the anapole moment is not a measurable, that is, 
not separately a gauge-invariant quantity. (However, we will later see that the 
dominate contribution to nuclear anapole moments is well defined and separately 
gauge invariant.) 
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Figure 1: Weak radiative corrections to electron-proton scattering include anapole 
contributions (a) as well as conributions that do not correspond to a single virtual 
photon exchange (b). 

Our primary focus in this review is the nuclear anapole moment contribution to 
atomic parity nonconservation (PNC). Exquisitely precise (sub 1%) measurements 
of atomic PNC have become possible in the past few years (Q). The primary focus of 
these studies has been to obtain precise values of the strength of direct Z° exchange 
between electrons and the nucleus. The PNC effects are dominated by the exchange 
involving an axial Z° coupling to the electrons and a vector coupling to the nucleus. 
The nuclear coupling is thus coherent, proportional to the nuclear weak vector 
charge (approximately the neutron number), and independent of the nuclear spin 
direction. It is widely recognized that these atomic measurements are important 
tests of the standard electroweak model and its possible extensions, complementing 
what has been learned at high energy accelerators that directly probe physics 
near the Z° pole §,|,@). The comparison between precision measurements at 
atomic energies and accelerator energies could reveal the subtle influence of new 
interactions beyond the standard model. 

At low energies one expects weak radiative corrections, including the anapole 
contribution, to interfere with the dominant tree-level weak amplitude, producing 
corrections to observables of relative size a ~ 1%, where a is the fine structure 
constant. Therefore identifying the anapole contribution in weak processes would 
seem a daunting task. The possibility that nuclear spin-dependent PNC measure- 
ments in heavy atoms might prove an exception was first pointed out two decades 
ago by Flambaum, Khriplovich, and collaborators @. Such spin-dependent PNC 
effects involve a vector coupling to the electrons and an axial coupling to the 
nucleus. In this case the contribution of Z° exchange is considerably reduced 
because the nuclear interaction is no longer coherent, but instead samples, in a 
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naive picture, only the spin of the last, unpaired nucleon within an odd-^4 nucleus. 
Additional suppression comes from the small Z° vector coupling to the electron, 
~ (4sin 2 #y^ — l)/2 ~ —0.05. In contrast, the anapole moment is enhanced in 
heavy nuclei, growing as A 2 / 3 , where A is the atomic number, and thus increasing 
in proportion to the nuclear surface area. The net result is the surprising con- 
clusion that the anapole contribution - the weak radiative "correction" - will in 
fact dominate over the tree-level direct Z° spin-dependent contribution for nuclei 
heavier than A ~ 20. 

The principal motivation for this review was the recent successful determina- 
tion of the nuclear spin-dependent contribution to atomic PNC in 133 Cs @. The 
extraction of this contribution from the much larger coherent PNC signal was 
accomplished by studying the dependence of the PNC signal on the choice of hy- 
perfine level. As the hyperfine differences are quite small, this extraction requires 
detailed control of possible sources of systematic error. Over a decade ago the Col- 
orado group ([h]) succeeded in measuring PNC effects in 133 Cs to 2.2%, providing 
a tentative identification of the anapole moment. Six years ago the Seattle (|iT|) 
and Oxford (12) groups achieved an accuracy of 1.2% and 2.9%, respectively, in 
PNC measurements in natural Tl, resulting in bounds on the anapole moment in 
this nucleus. Finally, in 1997 a definite measurement was made: the Colorado 
group of Wood et al. (||) reported PNC measurements in 133 Cs at the 0.35% level, 
from which a 7a nuclear spin-dependent PNC signal was extracted. This spin 
dependence is clearly well above the expected signal from Z° exchange alone but 
compatible with the enhancement expected from anapole effects. 

As we will describe in this review, a nucleus generates an anapole moment 
through the weak NN interaction, which operates at long range in the nucleus 
by meson exchange, where one meson-nucleon vertex is strong and the other 
weak (13, 14). Such interactions admix into the nuclear ground state odd-parity 
amplitudes and also induce new PNC nuclear currents. Thus the exquisite precision 
now achieved in atomic PNC studies has opened up a new window on the hadronic 
weak interaction. This interaction has proven more elusive than the weak interac- 
tions involving leptons. Whereas the charge-changing hadronic weak interactions 
can be studied in strangeness- or charm-changing decays, the standard model pre- 
dicts that neutral-current interactions do not change flavor. Thus hadronic inter- 
actions mediated by the Z° can only be studied in NN interactions, where PNC 
must be exploited to separate this contribution from much larger strong and elec- 
tromagnetic effects. This is a difficult task, and only a few hadronic experiments 
have achieved the requisite precision and even fewer (notably those done with p+p 
scattering or in certain light nuclei with special properties) can be interpreted, rea- 
sonably free of nuclear structure uncertainties. The importance of atomic anapole 
moment measurements is not only that they supplement the hadronic data, but 
also that they are sensitive to long-range pion-exchange PNC interactions where 
the effects of the neutral current can be isolated. We will see, however, that the ex- 
traction of weak coupling constant constraints from anapole moments are requires 
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some nonrivial nuclear physics analysis. 

The plan of this review is as follows. In Section 2 we discuss general properties 
of anapole moments, the associated current distribution, and the spin-dependent 
interaction nuclear anapole moments generate in atoms. In Sections 3 and 4 we 
summarize the experimental status of atomic anapole moment measurements and 
the progress theory has made in estimating the size of nuclear anapole moments 
from an underlying model of the hadronic weak interaction. In Section 5 we com- 
pare weak meson- nucleon coupling constraints obtained from the 133 Cs anapole 
result and Tl limits with similar constraints obtained from PNC nuclear observ- 
ables such as p + p scattering and the circular polarization of gamma rays emitted 
from 18 F. It becomes clear that, while all measurements are consistent with the 
broad "reasonable ranges" for PNC meson-nucleon couplings that theorists have 
defined (|l3|), there is some level of disagreement between various experiments when 
a global fit is performed. We conclude by discussing prospects for improving our 
experimental and theoretical knowledge of anapole moments. 

2 ANAPOLE MOMENTS AND ATOMS 

In this section we discuss general properties of nuclear anapole moments and the 
interactions they induce in atoms like 133 Cs and Tl. 

2.1 Electromagnetic Moments 

A standard multipole decomposition ( |i~5| ) of the electromagnetic current groups 
interactions according to their multipolarity and symmetry properties into charge 
(Cj), transverse electric (Tj l ), and transverse magnetic (T™ a9 ) operators, where 
J denotes the multipole rank. Static moments correspond to diagonal matrix 
elements of these multipole operators. From the transformation properties of the 
ordinary electromagnetic current operator under parity (P) and time-reversal (T) 

PJH2,t)p- x = Jp&t) 
fj»(x,t)f- l = J^x,-t), (1) 

and the constraint of hermiticity, it is readily verified that the non-zero moments 
arise from matrix elements of the even-J projections of Cj and the odd-J projec- 
tions of T™ 9 (0). 

More possibilities arise if one turns on the weak interaction, which introduces 
both parity- and CP/T- violating terms. Thus it is interesting to classify possible 
static moments more generally, without the assumption that the underlying Hamil- 
tonian conserves electromagnetic symmetries. The results in Table [l] show that, in 
addition to the ordinary monopole charge and magnetic dipole moments, two new 
dipole moments then arise. One of these, corresponding to the expectation value 
of the C\ multipole, requires both parity and CP/T violation. This is the electric 
dipole moment which, while allowed in the standard model due to CP violation in 
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Table 1: Transformation properties of electromagnetic moments 



J 


Cj 
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rpmag 
J 





P-even, T-even 
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P-odd, T-odd 


P- 


•odd, T-even 


P-even, T-even 


2 


P-even, T-even 


P- 


■even, T-odd 


P-odd, T-odd 


3 


P-odd, T-odd 


P- 


■odd, T-even 


P-even, T-odd 



the quark mass matrix and in the 9 term, has yet to be detected experimentally. 
The second, an PI moment, requires only parity violation. This is the anapole 
moment. 

These monopole and dipole moments must therefore arise in the most general 
expression for the matrix element of a conserved four-current for a spin-^ particle 

- <rV)7 5 - i^p-^wsMp). (2) 

The two vector terms define the charge Fi(q 2 ) and magnetic P2(<? 2 ) form factors. 
The axial terms that follow are the anapole and electric dipole terms, respectively. 
The anapole term reduces in the nonrelativistic limit to 

-^U<7 M -<?V)75 _> ^.f{a-qq-a) 

a(q 2 ) 
M 2 

showing that the current is transverse and spin dependent. We define the anapole 
operator as 

A 1X = a(0)a lx . (4) 

Weak radiative corrections that modify physical processes involving elementary 
fermions thus include anapole moment contributions. One interesting example 
is the Majorana neutrino, where the identity under particle-antiparticle conjuga- 
tion requires the magnetic and electric dipole moments to vanish, but permits an 
anapole moment (|I6|). Yet the interest in anapole moments would likely have re- 
mained largely theoretical were it not for the realization that they might play a 
significant role in atomic PNC expriments. 

2.2 Atoms and the Generalized Siegert's Theorem 

A generalization of the anapole operator is helpful for the case of a composite 
system like a nucleus, where the current operator and wave functions are modified 



q 2 a ±: (3) 
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by the interactions among the constituents, inducing parity admixtures in the 
ground state and multi-body currents. In the standard multipole decomposition 
the El projection of the current is expressed in terms of the multipole operator 

J ix (q) E -X -iVtoti? x (q), (5) 

where 

Tf x (l) = ~ I d 3 xf(x) • V x (h(qx)Y llx (n x )). (6) 
q J 

It is well known that current conservation places constraints on the matrix ele- 
ments of Tj l x (q). A familiar example is the long- wavelength limit of Tf x (q) gener- 
ated from the ordinary, parity-conserving electromagnetic current. The operator 
then is p/M, which is of order v/c, where v is the nucleon velocity. It can be 
shown that the exchange-current contributions to the vector three-current are also 
of order v/c. For realistic models of the nucleus which account for the interactions 
among the nucleons, in general we lack a prescription for constructing interactions 
and currents consistently - and for renormalizing them appropriately to take into 
account the limited Hilbert spaces employed in nuclear models. As a result, there 
will be errors in evaluations of Tf x , owing to the imperfect construction of the 
current, that are necessarily of leading order in the velocity, v/c. 

Siegert (|l^) showed that the situation could be greatly improved by exploiting 
the continuity equation 

V- J(x) = -i[H,p{x)\ (7) 

to rewrite Tf, in the long wavelength limit, entirely in terms of the charge operator. 
This generates the familiar dipole form of Tj l , proportional to wf, where u> is the 
energy transfer. The importance of this rewriting is that the charge operator, 
which is of order (v/c)°, has exchange current corrections only of order (v/c) 2 , or 
of relative size ~ 1%. Thus the Siegert's form of Tf, in which the constraints of 
current conservation are fully exploited, is a far more controlled operator for use 
in nuclear calculations. 

The analogous situation arises for the anapole moment. While there are a vari- 
ety of forms of the anapole operator that are equivalent up to terms that vanish 
by current conservation, these forms are not equivalent operationally in realis- 
tic calculations because of model violations of current conservation. The simple, 
long-wavelength form of Siegert's theorm doesn't address the question of moments 
because it generates an operator proportional to u, which then vanishes for a di- 
agonal matrix element. Fortunately the generalization of Siegert's theorem ( |l8|) 
exists: at arbitrary q one can write Tf(q) = Sj(q) + Rj(q), where all components 
of the electromagnetic current that are constrained by current conservation have 
been isolated in Sj and expressed as a commutator of the charge operator with 
the nuclear Hamiltonian. All such terms then vanish for a static moment. The 
resulting generalized Siegert's form oiTf appropriate for diagonal matrix elements 
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is 

ff[(q) f =X° --^= [ d 5 rr 2 [J lx (f) + V^(Y 2 (n r ) Ji(r)) 1A ]. (8) 

diagonal 9 y D7T J 

From Eqs. (5) and (8) we have the appropriate threshold form of the current 
operator. The anapole operator is defined relative to the current operator as in 
Eqs. (3) and (4), yielding immediately (|l9|) 

M2 /" ,3 2r 



A 1X = -— J d 6 rr z [Ji\^) + V2vr(y 2 (O r ) ® Ji(r))u] (9) 
Other forms of the anapole operator are more commonly used, for example (p0|) 

in = -7T / d 3 rr 2 J 1A (r). (10) 



Apart from the normalization, which is a matter of convention, this form will be 
equivalent to Eq. (9) only if the nuclear model is sufficiently simple that exact 
current operators can be constructed. This would be the case for a simple central 
potential model, such as an independent-particle harmonic oscillator, for which 
the appropriate current operator is that of a free particle. This would also be the 
case for an independent-particle model with spin-orbit and orbital potentials of 
the form a ■ t or i ■ £, provided that the additional contributions to the current 
obtained by minimal substitution, p — > p — eA into £ = fx p, are included in the 
calculation (^). But in other, more realistic treatments of the nuclear physics, Eq. 
(9) is the unique form that fully enforces the constraints of current conservation, 
regardless of the complexity of the current operators. 

Figure 2 gives a classical picture of the anapole moment as a current winding 
within the nucleus. Although the currents on the inner and outer sides of the torus 
oppose one another, there is a net contribution to Eq.(9) because of r 2 weighting of 
the current, leading to an anapole moment that points upward. (Note the sign in 
Eq. (9).) The current distribution drawn in Fig. 2 is odd under reversal of parity, 
as is the ordinary Jf™. Thus it is easy to see that the corresponding anapole 
moment is a parity-odd operator. If, however, the current has a chirality - a small 
"pitch" corresponding to a left- or right-handed winding that would signal PNC - 
a parity-even contribution to the operator would be induced. This is the analog of 
evaluating Eq. (9) with the axial-vector current induced by the weak interaction, 
leading to an operator component that would have a nonzero expectation value for a 
parity-conserving ground state. Similarly, the anapole moment associated with J^ 1 
would have a nonzero expectation value if the nuclear ground state contains parity- 
odd components induced by the weak PNC component of the AW interaction. 



2.3 The Spin- Dependent Atomic Interaction 



The potential felt by atomic electrons due to the nuclear anapole moment is gen- 
erated in the standard way by the dipole current - dipole current interaction. As 
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a 




Figure 2: A toroidal current winding that would correspond to a nonzero anapole 
operator in Eq. (9). 



the (f appearing in the nuclear current cancels the photon propagator for a static 
interaction, this weak interaction has the contact form 

G F 



w 



—=kol ■ loir] 
y/2 



(11) 



where I and p(r) are the nuclear spin and density, and where a is the Dirac 
matrix acting on the electrons. The superscript s.d. denotes that this is the nuclear 
spin-dependent contribution to the atomic PNC interaction. The portion of this 
interaction generated by the nuclear anapole moment is then 



^anapole 



AtraV2< > /e 



M 2 Gf 



< Hl/ll/ > 



(12) 



where 



denotes a matrix element reduced in angular momentum and a is the fine 

1)(2/ + 1). [Note 



structure constant. The reduced matrix element of / is 
that elsewhere in the literature another definition of k is more commonly used, 
one that associates the nuclear ground state with a single-particle level of orbital 
angular momentum i and spin /. For example, the k used in Ref. (|8|) is obtained 
by dividing ours by the factor (-l) /+1 / 2+£ (/ + l/2)/(/(/ + 1)).] 

We have already noted that K anapo i e is just one of the contributions to n: isolating 
this anapole contribution would be clearly easier if n anapo i e were the dominant such 
contribution. It was the important observation of Flambaum and Khriplovich (|^) 
that the A 2 ^ growth of nuclear anapole moments, already implicit in the r 2 weight- 
ing of the current in Eq. (9), would lead to such dominance in heavy atoms. It is 
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also clear that, if n anapo i e could be extracted from atomic measurements, relating 
this quantity to the underlying sources of hadronic PNC involves specifying the 
PNC currents that contribute to Eq. (9), as well as the PNC admixtures in the 
nuclear ground state wave function through which the ordinary electromagnetic 
current operator will generate nonvanishing matrix elements of A\\. We will defer 
this question to Section 4, focusing now on the experimental progress that has 
produced values and limits for k. 



3 EXPERIMENTS 



The idea of the anapole moment languished for many years after the early work of 
Zeldovich and Vaks. In 1973 Henley, Huffman, and Yu (22) - who were unaware 
of the Zeldovich paper - noted that the anapole moment would contribute to PNC 
observables in high energy electron scattering off the proton. Then a revival of 
interest in the anapole moment occurred in 1980 when its possible relevance to 
ongoing experiments on atomic PNC was noticed @. In these experiments the 
tiny parity- violating mixing of S and P states induced by the coherent, nuclear- 
spin-independent Z° coupling to the nucleus was being measured (23). As both the 
coherence and the increased electronic overlap with the nucleus lead to larger PNC 
effects in heavy atoms, the experiments focused on atoms with A > 100. The A 2 / 3 
enhancement of the anapole moment makes the spin-dependent contribution to 
atomic PNC dominant over other vector(electron)-axial(nucleus) sources of PNC in 
such heavy atoms. Spin-dependent effects of the expected size - a few percent of the 
spin-independent signal - could, in principle, be extracted from the experiments by 
precisely comparing the amount of mixing for two different hyperfine components, 
i.e., electronic transitions that differ only in the orientation of the nuclear spin. 

Nevertheless, the spin-dependent PNC effects are painfully small, corresponding 
to state mixings on the order of parts in 10 13 . That level of precision was even- 
tually reached in a cesium PNC experiment and nearly reached in experiments 
with thallium. The isolation of a definite spin-dependent contribution in cesium 
provided the first confirmation that anapole moments exist. 



3.1 The Cs Experiment: Techniques 

The experiment in cesium is described in detail in Refs. (p4|). The central idea is 
to exploit the highly forbidden 6S-to-7S transition in atomic cesium. The strength 
of the transition depends on the handedness within the excitation region. That 
handedness is defined by various applied electric, magnetic, and laser fields and 
reverses with appropriate reversals of those applied fields. This reversal forces 
any PNC component to the excitation rate to change sign, thereby altering the 
excitation rate for the transition. 

As the fractional change in the excitation rate is very small, a great deal of 
work must be done to achieve a signal-to- noise ratio sufficient to see any effect. 
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Figure 3: Partial cesium energy-level diagram including the splitting of the 5 states 
by the magnetic field. The case of 540-nm light exciting the (F = 3, = 3) level 
is shown. Diode lasers 1 and 2 optically pump all of the atoms into the (3,3) level, 
and laser 3 drives the 65^=4 (i^ e t) to 6Pf=5 transition to detect the 75 excitation. 
PNC is also measured for excitation from the (3,-3), (4,4), and (4,-4) 65 levels. 
The diode lasers excite different transitions for the latter two cases. 



Then even more work must be invested to verify that the detected effect is truly 
a violation of parity and not a spurious signal arising from systematic errors such 
as imperfect reversals or alignments of the fields that define the handedness of the 
experiment. 

In the absence of electric fields or parity violating interactions, the electric dipole 
(El) transition between the 65 and 75 states of the cesium atom (Fig. 3) is 
forbidden. As the nuclear spin of 133 Cs is I = 7/2, these S 1 / 2 levels combine with 
the nuclear ground state to form hyperfine states of total angular momentum F = 
4 and 3. A PNC interaction mixes a small amount (~ 1CP 11 ) of the neighboring 
6P3/2 and 7P3/2 states into the 65 and 75 states: the P 3 / 2 hyperfine levels have 
F = 2-5, so that PNC mixing with both the F = 3 and 4 hyperfine 5 states 
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takes place. The induced 75 <-> 6S 1 PNC -El amplitude is thus proportional to 
the product of the PNC mixing and El matrix elements coupling the S and P 
hyperfine levels (and inversely proportional to the energy differences between the 
S and P levels). To obtain a measurable observable that is first order in this 
small amplitude, a d.c. electric field E is applied that also mixes S and P states. 
This field generates an interfering "Stark-induced" EX transition amplitude Ap 
that is typically 10 5 times larger than the PNC-induced El amplitude Ap^c- A 
complete analysis of the relevant transition rates between the various hyperfine 



magnetic states (F,mp) is given in Ref. (24). These rates involve a straightfoward 



evaluation of the hyperfine matrix elements for the Stark amplitude, the nuclear- 
spin-dependent PNC Hamiltonian given in Eq. (11), and the spin-independent 
PNC Hamiltonian for the much stronger coherent interaction with the nucleus, 

= ^Qwlspir). (13) 

Here 75 is the axial coupling to the electrons, while the tree-level standard-model 
result for the vector weak charge of a point nucleus is 

Q w = Z(l - 4 sin 2 9 W ) ~ N ~ —N. (14) 

We omit the somewhat tedious angular momentum algebra here. To generate a 
nonzero interference between Ap^c and Ap, which differ by a relative phase of i, 
the 6S-to-7S transition must be excited with an elliptically polarized laser field of 
the form e z z + p i Im^) x, where the handedness of the polarization p is ± 1, e x is 
the component of the oscillating electric field parallel to the d.c. Stark field E, and 
e z is the oscillating field component in the direction perpendicular to E. The PNC 
contribution to the transition rate is measured on both the QSp^-to-TSp^ and 
6S'i^ = 4-to-75'p = 3 transitions. To a good approximation the only difference between 
the two transitions is the reversal of the nuclear spin. Thus simply taking the 
difference between the PNC contributions to the hyperfine transition rates isolates 
the nuclear-spin-dependent PNC coupling k of Eq. (11). 

For both transitions the atoms are initially populated and excited only out of 
states with extreme values of mp, (F = 3,mp = ±3) and (F = 4, mp = ±4). For 
these cases and the configuration of electric and magnetic fields shown in Fig. 4, 
the transition rate is 

R = \A E + A PNC \ 2 ~ ^ElelC^F.mp-F'.m'p) 

+ 4/3^e :c pIm(e x )Im(Elp NC )C 2 (F,m F ;F / ,mt) (15) 

where (3 is the tensor transition polar izability, E x = \E\ is the d.c. electric field 
strength, and C\ and Ci are combinations of Clebsch-Gordan coefficients whose 
detailed form we have suppressed (|24|). They depend on the initial- and final-state 
hyperfine labels {F,mp) and (F',m' F ) and transform as Ci(m^) — > Ci(-mp) and 
C2(mp) — > —C2(—mp) under reversal of the magnetic labels. The tiny contribution 



Anapole Moments 



13 



proportional to Ap NC has been neglected as well as the small 6S — 7S parity- 
conserving magnetic dipole transition amplitude Ami- (As discussed later, Ami 
can be a source of many systematic errors that must be addressed carefully.) 




HFP,ZP 

diode 
lasers 



laser 



Figure 4: Schematic of the cesium apparatus. In the interaction region B is along 
the z axis, E is along the x axis, and the 540-nm dye laser beam defines the y axis. 



(16) 



The experimental quantity of interest is the fractional PNC modulation in the 
transition rate 

AR _ 2Im(e x )Im(El PNC ) 
~R~ ~ e z (3E 

that modulates with the reversals of E, m, and p. In the experiment, there are 
actually five "parity" reversals because mi? is reversed in three different ways. 
The use of this large number (five) of independent reversals is essential for the 
detection and elimination of systematic errors, providing a great deal of redundancy 
for the PNC signal. This redundancy means that although no single reversal is 
perfect, the product of all the imperfections is far smaller than the PNC signal. 
Furthermore, the signal modulations accompanying various combinations of field 
reversals provide additional information about the field reversal imperfections and 
orientations that help identify potential systematic errors. 
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3.2 The Cs Experiment: The Apparatus 

A simplified schematic of the apparatus used by Wood et al. @ is shown in Fig. 
4. An effusive beam of atomic cesium is produced by a heated oven and collimated 
with a multichannel capillary array nozzle. The beam is optically pumped into 
the desired {F,mp) state by light from diode lasers 1 and 2 (24,p5|). The 2- 



cm-wide beam of polarized atoms intersects the 540-nm standing wave (Gaussian 
diameter, 0.8 mm) that is inside a high-finesse (100,000) Fabry-Perot power buildup 
cavity (PBC). The PBC not only enhances the transition rate, but the standing 
wave geometry also greatly suppresses the troublesome modulation arising from 
Ae\ — Ami interference (p4|). The 540-nm light originates from a dye laser whose 
frequency is tightly locked to the resonant frequency of the PBC by a high-speed 
servosystem (p6|). Before entering the PBC, the dye laser light passes through an 
intensity stabilizer, an optical isolator, and a polarization control system made 
up of a half-wave plate, Pockels cell, and adjustable birefringence compensator 
plate. The ellipticity of the light is controlled by rotating the half-wave plate, 
and the handedness of the ellipse is reversed by switching the sign of the 1/4 
(quarter wavelength) voltage that is applied to the Pockels cell. The intensity 
stabilizer holds the amount of light transmitted through the cavity constant and 
hence stabilizes the field inside of the cavity. This field corresponds to about 2.5 
kW of circulating power. The PBC resonant frequency is held at the frequency of 
the desired atomic transition by a servosystem that translates the input mirror. 

The d.c. electric field in the interaction region is produced by applying ~ 500 
volts between two parallel 5-cm by 9- cm conducting plates, separated by about 
1 cm. The plates are made of flat pieces of Pyrex glass coated with 100 nm of 
molybdenum. Both field plates were divided into five electrically separate segments. 
This division makes it possible to apply small uniform and gradient electric fields 
along the y-axis for auxiliary diagnostic experiments. The entire buildup cavity 
and field-plate mounting system is rather elaborate to ensure precise alignment 
and extreme mechanical stability. 

After being excited out of the populated 65 hyperfine level up to the 75 state, 
an atom will decay by way of the 6P states to the previously empty 65 hyperfine 
level (Fd e t) more than 60% of the time. The population of F^t is detected 10 
cm downstream of the interaction region using laser florescence. Light from diode 
laser 3 excites each atom in F^et to the 6P3/2 state many times. The resulting 
scattered photons are detected by a silicon photodiode that sits just below the 
atomic beam. When the 65^=3-to-75p=4 line is measured, the detection laser 
drives the 65i?=4-to-6i-3/2,.F=5 cycling transition (Fig. 3). About 240 photons 
per Fdet atom are detected. For the 65^=4-75^=3 line, the detection transition 
is the 65i?=3-to-6P 3 / 2) ir=2 cycling transition. This cycle gives about 100 detected 
photons per F^et atom. The signal-to-noise ratio for this transition is about 20% 
lower. During each half cycle of the most rapid field reversal (E) and after the 
switching transient has passed, the detector photocurrent is integrated, digitized, 
and stored. For each stored value, the computer also records the field and spin 
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orientations. 

The signal-to-noise ratio needed for this experiment puts extreme requirements 
on laser stability. A fluctuation in the intensity, frequency, or direction of the light 
from any of the four lasers will introduce noise in the detected atomic fluorescence. 
The most extensive control is needed for the dye laser ( |26| ) , but the requirements on 
the three diode lasers used for optical pumping and detection are also severe. These 
requirements motivated substantial development of diode-laser stabilization tech- 
nology (p7|). Both optical and electronic feedback were used to lock the frequency 
of each diode laser to the desired atomic transition using saturated absorption 
spectrometers. 

Extensive and precise control of magnetic fields was required in the experiment. 
In the optical pumping region, there is a uniform 2.5-G field that must point in 
the y direction (parallel to the pumping laser beams). In the interaction region, 
a 6.4-G field must point precisely along either ±z. Between the two regions, the 
magnetic field must rotate gently enough that the atomic spins follow it adiabati- 
cally. Finally, the field must be near zero in the detection region, and it is necessary 
to precisely reverse the fields in the optical pumping and interaction regions inde- 
pendently without significantly perturbing the fields in the other two regions. The 
setup required the use of 23 magnetic field coils of various shapes to provide the 
necessary fields and gradients. Most of these coils were driven with both reversing 
and nonreversing components of current. 

Many additional elements were required to achieve sufficiently precise alignment 
and control of all aspects of the apparatus. These include 31 different servosystems 
to ensure optical, mechanical, and thermal stability. 

3.3 The Cs Experiment: Data and Results 

Approximately 350 hours of data on the PNC observables were acquired in five runs 
distributed over an 8-month period. Each of the five runs followed the same basic 
procedure. First, a set of auxiliary experiments was carried out to measure and set 
numerous quantities: (i) all three components of the average E and B fields and 
their y gradients in the interaction region; (ii) the magnitude and orientation of the 
birefringence of the PBC output-mirror coating; (iii) the polarization-dependent 
power modulation of the green laser light; and (iv) the populations of the mp levels 
of the atomic beam as it entered the interaction region. After these measurements 
were completed, the four laser frequencies were locked to the desired hyperfine 
transitions and data were acquired in blocks of about 1.5 hours each. During this 
time, five parity reversals - the electric field, laser polarization, and three ways 
of reversing the mp state being excited (reversing the polarization of the optical 
pumping light of laser 2, reversing the optical pumping magnetic field relative to 
the pumping light, and reversing the magnetic field in the 6S-7S excitation region) 
- were carried out at different rates. The electric field was reversed at 27 Hz and 
the others at various lesser rates. The relative phases of the various reversals were 
regularly shifted by one half cycle. Before and after each of the 1.5-hour blocks, the 
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polarization of the 540-nm standing-wave field was measured and set. At regular 
intervals, the PBC output mirror was rotated by 90.0 degrees, and at irregular 
intervals, the frequencies of the four lasers were changed to measure PNC on the 
other 65- 75 hyperfine line. At the end of each data run (20 to 30 1.5-hour blocks of 
data), the initial auxiliary experiments were repeated to check that the quantities 
described above had not changed significantly. 

The typical size of the 65 — 75 signal from the photodiode was 200 nA on the 
3-4 line and 85 nA on the 4-3 line. These measurements corresponded to about 
0.5% of the atomic beam undergoing a 65 — 75 transition. The signal-to-noise 
ratio for measuring the 65 — 75 transition rate was typically 55,000/V Hz, and 
the parity violating modulation was typically about six parts in 10 6 . Because the 
technical noise was small and many photons were detected from each atom that 
had undergone a 65 — 75 excitation, the noise was dominated by the shot-noise 
fluctuations in the number of atoms making the 65 — 75 transition. Data was 
taken over a range of polarization ellipticities and electric fields. 

The data analysis used to find the PNC modulation for each block of data was 
relatively simple. Appropriate combinations of fractional differences in the signal 
sizes were calculated for each of the 32 different field configurations to find the 
fraction of the rate that modulated with all five reversals. From the measured 
values of electric field and laser polarization this fractional modulation was then 
converted into Im(ElpNc)/P- Various other small calibration corrections were 
required. 

The major concern in this experiment was possible systematic errors arising 
from spurious signals that modulate under all five parity reversals, thus mimicking 
PNC. Roughly 20 times more data were taken in the investigation and elimination 
of such errors than in the actual PNC measurement. Several small errors associated 
with stray and misaligned fields were encountered, as in previous PNC measure- 
ments (24), and were treated as before. An exhaustive analysis was carried out of 
all possible combinations of static and oscillating electric and magnetic fields that 
could mimic a PNC signal. All of the stray (defined as nonreversing) and mis- 
aligned d.c. electric and magnetic fields and their gradients, and many of the laser 
field components, could be determined by looking at appropriate modulations in 
the 65 — 75 rate under various conditions. Many of these quantities were extracted 
from the 31 different modulation combinations observed in real time while taking 
PNC data, and the remaining components were determined by the auxiliary exper- 
iments that were interspersed with the PNC runs. Many tests were also performed 
to ensure the necessary stability of the relevant fields. 

Although the procedures were similar in concept to previous work, the Wood 
et al. (||) measurements were more difficult and time consuming due to the higher 
accuracy required. It proved necessary to consider not only the average fields, 
but also their gradients across the interaction region. Several small spurious sig- 
nals were identified and removed. The absence of any systematic effects from the 
troublesome Ae — Ami interference was confirmed by the independence of results 
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for various laser polarizations sensitive to this interference. Other cross-checks 
included enhancing sources of error to confirm the predicted response and per- 
forming data analyses to verify that transition rate variations were consistent with 
fundamental shot-noise fluctuations. 
The final results are 

- Im(El PNC )//3 = 1.6349 ± 0.0080 mV/cm 6S F=4 <-» 7S F=3 

= 1.5576 ± 0.0077 mV/cm 6S F=3 <-» 7S F=4 (17) 

yielding for the nuclear-spin-dependent difference of interest for the anapole mo- 
ment 

Im(Elp^- c )//3 = 0.077 ± 0.011 mV/cm, (18) 

a 7<T effect. The statistical uncertainties for the two transitions, 0.0078 and 0.0073 
mV/cm, respectively, dominate the error. The systematic uncertainties are based 
on statistical uncertainties in the determination of various calibration factors and 
systematic shifts, and therefore, it is appropriate to add them in quadrature. The 
final results are in good agreement with previous measurements in cesium but are 



much more precise. From the best available atomic calculation (28, 29) one then 
finds 

k( 133 Cs) = 0.112 ±0.016, (19) 
a result we will find is dominated by the anapole moment. 



3.4 The Tl Experiment 

Experiments measuring parity violation in atomic thallium (70.5% 205 T1, 29.5% 
203 T1) have not yet detected a nuclear spin dependent/anapole moment contribu- 
tion, but they have achieved accuracies very near the level where anapole effects are 
expected. We will see that the resulting limits are interesting from the perspective 
of hadronic PNC. An interference of parity- allowed and PNC contributions to an 
atomic transition is observed, just as in cesium. However, in thallium the transition 
is an allowed magnetic dipole transition, 6P1/2 ^ 6P 3 /2- The PNC interference is 
between the allowed Ml transition amplitude and the PNC El amplitude arising 
from weak interaction effects that mix S states into the P1/2 ground state. As Tl 
has an / = 1/2 nuclear ground state, the P\/2 state has F = and F = 1 hyperfine 
sublevels which will be affected differently by nuclear-spin-dependent sources of 
PNC. Thus the search for spin-dependent PNC effects requires a comparison of 
the strength of PNC transitions out of the F = and F = 1 states. Groups at the 
University of Washington (|Tl"| ) and Oxford (|i~2] ) have made such comparisons. 

The optical rotation of linearly polarized light that is produced by the interfer- 
ence of El and Ml transitions is measured in a large vapor cell of atomic thallium, 
as shown schematically in Fig. 5. A beam of linearly polarized laser light passes 
through the vapor and then through a nearly crossed linear polarizer followed by a 
sensitive detector. A Faraday modulator rotates the plane of polarization back and 
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vapor cell detector 




empty cell 



Figure 5: Schematic of the thallium vapor cell apparatus. 



forth through the perfectly crossed position. This allows phase-sensitive detection 
that minimizes sensitivity to drifts in the signal baseline. The laser is scanned 
repeatedly across the transition and the detected light signal is then fit to the pre- 
dicted lineshape. The PNC rotation reverses sign as the laser is tuned through line 
center. From the fit to the dispersion-like line shape the amplitude of the rotation 
can be determined. To distinguish the signal of interest from spurious frequency- 
dependent rotations in optics, the vapor cell is interchanged with a "dummy" cell 
that has identical optics but no atomic vapor. An example of the data obtained in 



this fashion is shown in Fig. 6, taken from Ref (11). 

It can be seen that the agreement between predicted and observed signals is 
excellent. As in the cesium experiment hundreds of hours of data were taken and 
considerable effort was made to eliminate possible systematic errors. The thallium 
oven was carefully constructed and shielded to avoid magnetic fields that would 
cause spurious Faraday rotation. Data were acquired over a wide range of pressures 
of thallium and a variety of tests of the calibration were carried out. 

The final results of the experiments are expressed as the ratio of parity vio- 
lating El to magnetic dipole Ml transition amplitudes, R = Im(ElpNc/Ml). 
The nuclear-spin-independent PNC signals obtained are R s l ' = (—14.68 ± 0.17) x 
1(T 8 (H) (normalized to 205 T1) and (-15.68 ±0.45) x 10~ 8 @, corresponding to 
1.2% and 2.9% accuracy. The nuclear-spin-dependent effects are consistent with 
zero, 

R sA -{T\) = (0.15 ± 0.20) x 10" 8 Seattle 

= (-0.04 ± 0.20) x 10~ 8 Oxford. (20) 
The resulting constraints on k are 

re(Tl) = 0.29 ±0.40 Seattle 
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Figure 6: Combined PNC optical rotation for the thallium data cycles (a total of ~ 
230 h) for F = 1 (points) compared to the theoretical (dispersion x transmission) 
line shape (solid line). 



= -0.08 ±0.40 Oxford. (21) 

The 133 Cs and Tl spin-dependent measurements involve odd-proton nuclei, which 
leads to the prediction that they test a similar combination of isospin components 
of the weak hadronic PNC potential, as we will see below. Note that an anapole 
measurement for an odd-neutron nucleus would test a roughly orthogonal combina- 
tion of isoscalar and isovector contributions to the hadronic PNC potential. Several 
atomic PNC efforts are underway that could produce such constraints (3C,|3l|). 



4 ANAPOLE MOMENTS AND HADRONIC PNC 



In this section we discuss the extraction of n an apoie from k, the relation between 
^anapole and the weak hadronic interaction, and the mechanisms by which a nucleus 
generates an anapole moment. 
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4-1 Extracting the Anapole Moment 

We have seen that atomic PNC measurements place the following constraints on 
the strength of nuclear-spin-dependent electron-nucleus interaction 

k( 133 Cs) = 0.112 ±0.016 

k(TI) = 0.29 ±0.40. (22) 

where we have used the Seattle Tl result above because, due to its central value, 
it proves to be the more restrictive in the weak meson-nucleon parameter region 
favored by nuclear PNC experiments. The three principal contributions to n 

K = ^anapole + K Z° + K Q W (^3) 

arise from the nuclear anapole moment, the vector(electron)-axial(nucleus) tree- 
level Z° exchange, and a term generated by the combined effects of the coherent 
Z° and magnetic hyperfine interactions between the electrons and the nucleus. In 
the theory discussions below we will treat the Tl constraint as one on the principal 
isotope, 205 T1 (70.5%). The other isotope, 203 T1 (29.5%), differs in structure only 
by a pair of neutrons, and thus should have very similar properties. 

The tree-level vector(electron)-axial(nucleus) Z° interaction generates a contri- 
bution 

K z° = — - 4 sin Ow) 7777— , (24) 

2 </||7||7> 

where g& = 1-26 is the axial- vector coupling and sin 2 9w = 0.223. To get a 
rough feel for this contribution, we can evaluate the nuclear matrix element in the 
extreme single-particle limit. 133 Cs would then be described as a unpaired 1*77/2 
proton outside a closed core, while 205 T1 corresponds to an unpaired 3s 1/2 proton. 
In this limit 

jingle particle = _ ( _ 1) /-^-m t _ 4 sin 2 9 W ) (25) 

where I is the single-particle orbital angular momentum and nit the z-component 
of isospin (with m t = 1/2 denoting a proton). Thus the single-particle estimates 
for 133 Cs and 205 T1 are 0.0151 and -0.136, respectively. Nuclear models of various 
types have been employed to try to estimate the effects of strong correlations in 
quenching the Gamow- Teller matrix element in Eq. (24) from its single- particle 



value (TL9, 32, 33, 34). We will employ the results of a recent large-basis shell model 



(SM) study (|!) here, which yielded 

4f( 133 Cs) = 0.0140 

4f( 205 Tl) = -0.127. (26) 

(We will describe these calculations in more detail later.) In addition, one-loop 
standard model electroweak radiative corrections, which we do not include in the 
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numerical results below, modify the tree-level expression in Eq. (24) somewhat, 
reducing the isovector contribution and inducing a small isoscalar amplitude (35). 



A second contribution to k is generated by the combined effects of the coherent 
vector Z° coupling to the nucleus, proportional to the weak charge Qw, and the 



magnetic hyperfine interaction (36). From the measured nuclear weak charge and 



magnetic moment Bouchiat and Piketty find (|3 

kq w { 133 Cs) = 0.0078 

«Qw( 2 ° 5 Tl) = 0.044. (27) 

Thus the experimental values for the anapole contributions to k are obtained by 
subtracting the results of Eqs. (26) and (27) from Eqs. (22), yielding 

KanapoU^Cs) = 0.090 ± 0.016 

Kanapole ( 2 ° 5 T1) = 0.376 ± 0.400. (28) 

4-2 The Hadronic Weak Interaction 

The various mechanisms operating within the nucleus to generate K ana p i e arise 
from the hadronic weak interaction. As we noted in the introduction, the only 
practical strategy for studying the effects of Z° exchange between quarks is the 
investigation of the PNC NN interaction. As anapole moments are now measure- 
able in high precision atomic PNC experiments, they now become a part of that 
strategy. 

The low-energy hadronic weak interaction can be described by a phenomenolog- 
ical current-current Lagrangian (|i~4| ) 

L= < ^={J ] w J w + J ] z J z ) + h.c. (29) 

where Jw and Jz are the charged and neutral weak currents, respectively. If one 
considers only the light-quark (u, d, s) contributions to these currents, then Jw has 
two components 

Jw = cos 0c Jw + sin OcJw (3°) 

where sin^c* ~ 0.22. The current J w drives the u — > d transition and transforms 
as AI = 1, AS" = 0, while J w drives the u — > s transition and transforms as 
AI = 1/2, AS = —1. (Here I denotes isospin and S strangeness.) The neutral 
current also has two components, and J^, which transform as AI = 0, AS = 
and AI = 1, AS = 0, respectively. The AS = NN interaction is then governed 
by the following piece of Eq. (29) 

l as=o = 9^ [cos *e c {j w )U w + W6 c (j w )U w 

+ W + (4) t 4 + (Jz) , 4 + W]H 6 (si) 
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An important aspect of this equation is its isospin content. The symmetric 
product of two Jyy (AI = 1) currents transform as AI = and 2, while the 
symmetric product of two (AI = 1/2) currents transforms as AI =1. Thus the 
AI = 1 component of the charged current weak NN interaction is suppressed by 
tan 2 9c compared to the AI = 0,2 components, while the neutral current AI = 1 
contribution is unsuppressed. It follows that the neutral current should dominate 
this isospin channel. 

The main goals of the field have been to isolate this neutral current and, more 
generally, to understand the mechanism by which the weak force is communicated 
over the relatively long distances characterizing low-energy NN scattering or NN 
interactions within the nucleus. While the standard model specifies the elementary 
couplings of the weak bosons to quarks, these vertices are dressed by the strong in- 
teraction to form couplings between physical particles, such as mesons to nucleons. 
We know, from the AI = 1/2 rule in strangeness-changing hadronic weak interac- 
tions (the strong enhancement of the ratio of AI = 1/2 to 3/2 amplitudes), that 
strong interaction effects can substantially alter couplings from their underlying 
bare values. 

The low-energy NN weak interaction is conventionally described in a one-meson- 
exchange model, where one meson-nucleon vertex is weak and the other strong: this 
long-distance mechanism dominates at nuclear densities. Six weak couplings, f w , 
hp, hp, h 2 p , h^, and hjj, characterize the strengths of the isovector tt, isoscalar/isovector/isotensor 
p, and isoscalar /isovector uj weak meson-nucleon couplings (|l3|). This model is not 
as restrictive as it may first appear. First, CP invariance forbids any coupling 
between neutral J = mesons and on-shell nucleons, eliminating a number of 
candidate interactions (^). Second, the most general low-energy PNC interaction 
contains only five independent S — P amplitudes. From this perspective, the de- 
scription in terms of ir, p, and uj exchange can be viewed as an effective theory, 
valid at momentum scales much below the inverse range of the vector mesons. At 
low momentum the detailed short-range behavior of the S — P potentials is not 
resolvable: thus one could characterize the short-range weak interaction by five 
contact interactions corresponding to these independent S — P amplitudes, sup- 
plemented by long-range it exchange. The six meson-nucleon weak couplings allow 
one to mimic these six degrees of freedom. 

Denoting the isoscalar strong meson-nucleon couplings by g n NN, 9 pi and g u , the 
resulting NN PNC potential is @ 



V PNC (r) = ^[f(l)xf(2)] 3 P(l) + ,l2)].u(rl 
g p fh p r(l)-r(2) + %rs(l)+r 3 (2)] 



1 

2M 



^[3T3(l)r 3 (2)-f(l).f(2)] 
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Table 2: Weak meson- nucleon coupling "best values" and "reasonable ranges" (in 
units of 10 -6 ) from the standard-model calculations of Desplanques, Donoghue, 
and Holstein 



Coupling "Best Value" "Reasonable Range 1 





0.46 


0.00- 


♦1.14 


s 


-1.14 


1.14- 


-3.08 


-0.019 


-0.038- 


^0.00 


-0.95 


-0.76- 


♦-1.10 


-0.19 


0.57- 


-1.03 


K 


-0.11 


-0.19- 


>-0.08 



x ((1 + » y )i[a(l) x <r(2)] ■ u p (r) + [<r(l) - a(2)\ ■ v p (r)) 

+ 9u>(hl + ^[T 3 (l)+T 3 (2)]\ 

x ((1 + Ha)i[ff(l) x a(2)} ■ u w (r) + [<?(1) - a(2)} ■ v w (r)) 

+ ~[r 3 (l) - r 3 (2)][a(l) + B{2)\ ■ \g^h l J u {r) - g p h l p v p (r)} 



(32) 



where r = f\ — r2, u = [p, e~ mr /4-Trr], v = {p, e~ mr /4irr}, and p = Pi - P2, with 
r*i and p\ the coordinate and momentum of nucleon 1. This expression is usually 
evaluated assuming vector dominance, which fixes the strong scalar and vector 
magnetic moments, fi s = -0.12 and \i v = 3.70. 

"Best values" and broad "reasonable ranges" for the weak meson-nucleon cou- 
plings were defined some time ago by Donoghue, Desplanques, and Holstein ( |i"3| ) 
(DDH), who deduced standard- model estimates for these vertices by using tech- 
niques like factorization, the quark model, and current algebra and sum rule meth- 
ods. The broad "reasonable ranges" reflect the large degree of uncertainty implicit 
in such approximate tools, as well as the potential consequences of missing physics, 
such as strange-quark amplitudes (|39|) . Nevertheless, the DDH results in Table ^ 
have provided experimentalists with benchmarks for PNC experiments. 

In an ideal world one would determine the low-energy NN S — P amplitudes, or 
equivalently the six weak meson-nucleon couplings, by a series of NN scattering ex- 
periments. Such experiments require measurements of asymmetries ~ 10 -7 , which 
is the natural scale for the ratio of weak and strong amplitudes, ^kGfvti^I g\ NN . 
Only a single iV./V measurement, the longitudinal analyzing power A z for p + p, 
has been successful (4C, P],P^). (Experiments have been done at 13.6, 45, and 221 



MeV.) These results have been supplemented by a number of PNC measurements 
in nuclear systems, where accidental degeneracies between pairs of opposite-parity 
states can produce, in some cases, large enhancements in the PNC signal. The 
experiments include A z for p + a at 46 MeV (]43|), the circular polarization P 7 of 
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the 7-ray emitted from the 1081 keV state in 18 F (fB|), and ^4 7 for the decay of 
the 110 keV state in polarized 19 F (fl5|). It is widely agreed that these experiments 
can be interpretted relatively free of nuclear structure uncertainties: the systems 
are either few-body, where quasi-exact structure calculations can be done, or in- 
volve special nuclei where the PNC mixing matrix elements can be calibrated from 
axial-charge (5 decay (|46|). An anaysis of these results, which have been in hand for 
some time, suggests that the isoscalar PNC NN interaction - which is dominated 
by p and uj exchange - is comparable to or slightly stronger than the DDH "best 
value," while the isovector interaction - dominated by tt exchange - is significantly 
weaker (< 1/3) (jlj). As the isovector channel is expected to be enhanced by neu- 
tral currents, there has been great interest in confirming this result. The 133 Cs 
anapole result thus provides a possible crosscheck on this tentative conclusion that 
an isospin anomaly, superficially like the AI = 1/2 rule in flavor-changing decays, 
may exist in the AS = weak NN interaction. 



4-3 Weak Meson-Nucleon Couplings and the Anapole Moment 

We denoted in the introduction that an anapole moment - unlike the magnetic 
moment - of elementary fermions is not a gauge invariant quantity. The very 
nice discussion of this point offered by Musolf (^TJ) ends with a very straightfor- 
ward explanation: because PNC corrections to the electromagnetic current couple 
only to virtual photons, the amplitude for PNC photon emission is not a physical 
amplitude and thus need not be gauge independent. However the long-distance 
contributions to the anapole moment in the nucleus - the meson cloud contribu- 
tions to the nucleon anapole moment and the many-body contributions due to the 
PNC NN interaction and associated exchange currents - are both the dominant 



contribution to the nuclear anapole moment and separately gauge invariant (19). 
These contributions, associated with the weak meson-nucleon couplings, will now 
be discussed. 

To evaluate the nuclear anapole moment in the context of some model of hadronic 
PNC one must take the ground-state expectation value of the operator in Eq. (9). 
If one works to first-order in the weak interaction, two types of terms contribute to 
the current matrix element. First are terms corresponding to the PNC weak contri- 
butions to J\\. In the context of the weak meson-nucleon couplings, these include 
both one-body terms - mesonic loop corrections to the ordinary electromagnetic 
current involving one weak and one strong vertex, as illustrated in Fig. 7a - and 
exchange currents (Fig. 7b), where the weak and strong vertices attach to differ- 
ent nucleons. We stress that if the anapole operator of Eq. (9) is employed, the 
exchange currents are model dependent in the sense that all constraints imposed 
by current conservation are already explicitly enforced. The argument sometimes 
heard that minimal substitution in simpler, independent-particle treatments some- 
how accounts for the exchange currents is incorrect: once Eq. (9) is employed, 
the surviving exchange currents depend on aspects of the NN interaction not con- 
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Figure 7: The one-body (a), exchange current (b), and nuclear polarization (c) 
contributions to the nuclear anapole moment. 



The second class of contributions comes from evaluating Eq. (9) for the ordinary, 
parity-conserving Jf™, which contributes through PNC admixtures in the nuclear 
wave function, as illustrated in Fig. 7c. The polarization term, which depends 
on the matrix elements of y PNC between the ground state and a complete set of 
excited opposite-parity states, can be enhanced if there exist opposite-parity states 
very near the ground state: the admixing is clearly inversely proportional to the 
energy difference (^7|). However, in the absence of such "accidental" degeneracies, 
one expects the mixing to be dominated by the giant resonances, the collective 
states at ~ 15 — 20 MeV in heavy nuclei that account for most of the El and other 
hrst-forbidden nuclear response. 

We discuss each of these contributions below, depending rather heavily on the 
recent work of Ref. (|3^). In Ref. (|32] ) an effort was made to model the various 
contributions to the anapole moment using the formalism that has become standard 
in other studies of hadronic PNC: the DDH weak-meson couplings, the two-body 
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yPNC Dasec [ on those couplings, large-basis SM nuclear wave functions, and a 
standard correlation function to modified the short-distance behavior of the SM 
two-nucleon density. Thus this treatment will allow us, in Section 5, to compare 
anapole, AW, and nuclear constraints on hadronic PNC on an equal footing. 

a) Nucleon anapole moment. The one-body PNC electromagnetic current derived 
from meson loop corrections (e.g., Fig. 7a) yields in the nonelativistic limit 

A 

A nucleon,c = £ ^ + ^(Q)^)] 0~\\(t). (33) 
i=l 

That is, this term is the sum of the anapole moments of the individual nucleons. 
As the contribution of spin-paired core nucleons cancel, one expects this term to 
depend on the unpaired valence nucleon in the odd-j4 nuclei of interest. Thus 
the one-body contribution should be roughly independent of A though - like the 
closely related example of magnetic moments - its value will depend on the shells 
occupied by the valence nucleon. Some time ago the pion contribution to a s (0) and 
a„(0) was evaluated (|l9|), yielding results proportional to ef^g^NN- The isoscalar 
coupling a s (0) proved to be about four times larger than a„(0). Recent extensions 
of this work have included the full set of one-loop contributions involving the DDH 
vector meson PNC couplings, using the framework of heavy baryon chiral pertur- 
bation theory and retaining contributions through 0(1/ A 2 ), where A x = 47ri ? 7r ~ 1 
GeV is the scale of chiral symmetry breaking (^). The contributions due to f n 
are consistent with the earlier work - the new a s (0) is about 1.3 times larger, while 
a v (0) is zero in this order - so that the dominance of a s (0) in Ref. ( |l9|) was is ex- 
plained as a consequence of heavy baryon chiral perturbation theory. The addition 
in Ref. ( |35| ) of the heavy meson contributions greatly enhances a v (0). An evalua- 
tion with DDH best value couplings yields a v (0) ~ 7a s (0). Thus the inclusion of 
heavy meson contributions substantially enhances the one-body anapole terms and 
alters the isospin character, generating opposite signs for the proton and neutron 
anapole moments. 

This then determines A™ ebody . The SM studies of Ref. @, which we will 
describe in more detail below, give matrix elements (J| | J2i=i 1-0 = -2-37 and 
2.53 and (J|| Y,f=i = -2.30 and 2.28 for 133 Cs and 205 T1, respectively. 

This leads to the nucleonic anapole contributions shown in Table [3|. 

b) Exchange currents. Two-body PNC currents arise from meson exchange dia- 
grams where the photon couples to the meson "in flight" (transition currents), or 
where the photon and either the PC or PNC meson-nucleon vertex creates and 
annihilates an AW pair (pair currents), as illustrated in Fig. 7b. The pionic ex- 
change currents, which have the longest range, were evaluated in Ref. ([H]). The 
explicit form of the pion contribution to Af^ cur is given there. Rather than quote 
the result, which is a bit complicated, we instead give the one-body reduction of 
the pair current contribution to that operator, which illustrates the underlying 
physics much more clearly. 

If one views a nucleus as a single-particle outside of a closed core, a two-body 
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Table 3: Shell model estimates of the anapole matrix element (I\ \A\\ \I)/e, ex- 
pressed as cofficients times the indicated weak couplings 



Nucleus 


Source 


u 


h u 

n P 


K 


ti 2 




hi 


i; «Cs 


nucleonic 


0.59 


0.87 


0.90 


0.36 


0.28 


0.29 




ex. cur. 


8.58 


0.02 


0.11 


0.06 


-0.57 


-0.57 




polariz. 


51.57 


-16.67 


-4.88 


-0.06 


-9.79 


-4.59 




total 


60.74 


-15.78 


-3.87 


0.36 


-10.09 


-4.87 


205 rpj 


nucleonic 


-0.63 


-0.86 


-0.96 


-0.35 


-0.29 


-0.29 




ex. cur. 


-3.54 


-0.01 


-0.06 


-0.03 


0.28 


0.28 




polariz. 


-13.86 


4.63 


1.34 


0.08 


2.77 


1.27 




total 


-18.03 


3.76 


0.33 


-0.30 


2.76 


1.26 



operator can be replaced by an equivalent one-body effective operator 

< a\A e / x f \f3 >= ]T < a5\A\l cur \(35 - 6(3 >, (34) 

8<F 

where the sum is taken over the nuclear core. Completing this sum in Fermi gas 
model, assuming a spin-symmetric but isospin-asymmetric core, yields the pionic 
pair current effective operator 

M A 

2.74a:(0)—Y / p(n)rt(* lx (i) + V^[Y 2 (r t )®a(i)} lx ) 

n 1=1 

(|u;f (1 " |ts(«)) + ^(1 + |ts(i))) , (35) 

where p(ri) is the nuclear density operator and cu^ (wjf) a proton (neutron) Fermi- 
gas response function that depends on kijkp, the nucleon momentum as a fraction 
of the Fermi momentum. The los vary only gently, ranging from 0.33 to 0.19 as 
ki/kF increases from to 1. Thus a suitable average value is ~ 0.25. 

The overall strength of A^(f' NN is given in terms of the pionic contribution to 
the single-particle anapole moment 

o j(0) ~ -1.6 ef * 9 * NN (36) 

to allow an easy comparison with the one-body current. Using a nuclear density 
of 0.195/fm 3 , one finds that the pair isoscalar ir exchange current then scales as 
~ 0.9yl 2 / 3 aJ(0). Unlike the nucleonic contribution, the exchange current contri- 
bution grows as A 2 / 3 . Clearly it will increasingly dominate over the nucleonic 
contribution as A increases. The isovector effective operator is smaller by a factor 
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of 2{Z — N)/3A, reflecting a cancellation between contributions from core protons 
and neutrons. 

The calculations of Ref. ( |l9| ) employed the full two-body form for the pionic 
currents, evaluating these from the SM two-body density matrices. A short-range 
correlation function was introduced to mimic the effects of hard-core correlations 
on the density matrix. The extension to include the p and u PNC couplings is a 
formidable task requiring evaluation of the p and u pair currents and the ppj and 
pirj currents. This calculation was only recently completed (|32[). The pp^f and 
pirj transition currents and the component of the oj pair current where the photon 
and PNC u> couplings are on different nucleon legs were found to be negligible, 
well below 1% of the dominant ir currents; the remaining heavy-meson currents 
are more important but, as can be seen in Table ||, still are suppressed relative to 
the pionic currents. The tabulated results were obtained from the SM calculations 
described below. 

c) Nuclear polarization contribution. The nuclear polarization contribution (Fig. 
7c) to the anapole moment is given by 

E (3^MME^R + h . c _ (37 ) 

where Af m is generated from inserting the ordinary electromagnetic current op- 
erator into Eq. (9), |/) is a ground state of good parity, H PNC is the PNC NN 
interaction, and the sum extends over a complete set of nuclear states n of an- 
gular momentum / and opposite parity. The extended Siegert's theorem again 



determines the form of Af 71 (19). 

Dmitriev, Khriplovich, and Telitsin ( |4g| ) have provided a series of estimates of 
the polarization contribution in the context of single particle models, using a one- 
body effective V PNC . The models include one of uniform density and ones based 
on the harmonic oscillator and Woods-Saxon potentials. The former yields an 
anapole operator that is explicitly proportional to the nuclear density - a quantity 
approximately constant in heavy nuclei because of nuclear saturation - and grows 
as A 2//3 , like the exchange current contribution. Numerically, however, the polar- 
ization contribution is about a factor of five larger than the exchange currents: it 
is the dominant contribution to axial PNC in heavy nuclei. 

For the present analysis we need a treatment of polarization contribution that 
begins with the DDH form of V PNC . This presents two nuclear structure chal- 
lenges. The first is the construction of a reasonable model for the ground state. 
The SM calculations of Refs. (HH) for 133 Cs were done in the canonical space 
between magic shells 50 and 82, lgj/2 — 2d 5 / 2 — — ^ s i/2 ~ 2^3/2- The protons 

were restricted to the first two of these shells and neutron holes to the last three, 
producing an m-scheme basis of about 200,000 . 205 T1 is described as a proton hole 
in the orbits immediate below the Z=82 closed shell (3si/2 — 2^3/2 — 2^5/2) coupled 
to two neutron holes in valence neutron space between magic numbers 126 and 82 
(3pi /2 - 2/ 5/2 - 3p 3/2 - H13/2 - 2/7/2 - 1^9/2)- 
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The second challenge is the completion of the intermediate state sum. Direct 
or moments methods ( |i~9| ) are impractical because of the dimensions of the spaces. 
However, because no nonzero El transitions exist among the valence orbitals, an 
alternative of completing the sum by closure, after replacing 1/E n by an average 
value (1/E), is quite attractive: the resulting product of the one-body operator 
Af 71 and two-body V PNC contracts to a two-body operator, so that only the two- 
body ground state density matrix is needed. The closure aproximation can be 
considered an identity if one knows the correct (1/E): in practical terms, this 
means demonstrating that it can be related reliably to some known quantity, with 
the position of the El giant resonance being the obvious candidate. A series of SM 
calculations has been completed for a series of light nuclei, where both the 1/E— 
and non-energy- weighted sums could be done (|32"|). It was found that there was a 
consistent relationship between the closure energy and the giant dipole position, 
though the three isospin contributions to y PNC must be treated separately. The 
average excitation energies found (as fractions of the dipole energy) are 0.604 ± 
0.056 (/#,/#), 0.899 ±0.090 (/„.), and 1.28 ± 0.14 (h 2 p ). The larger (1/E) for h° p 
and enhances these contributions to the anapole polarizability. The assumption 
that these relations also hold in the heavier nuclei 133 Cs and 205 Tl then fixes 
the appropriate average excitation energies for these nuclei. From the SM two body 
density matrix, the contracted two-body effective operator that results from the 
closure approximation (this includes the effects of a short-range correlation function 
on the two-nucleon matrix elements of V PNC ), and these excitation energies, one 
obtains the SM polarization results in Table ||[ 

Note that there are situations - chance ground-state parity doublets or nuclear 
octupole deformation - that can lead to enhanced polarizabilities. While such 
enhancement is unexpected in 133 Cs as the nucleus is approximately spherical, en- 
hancements have been discussed in connection with other nuclei of interest because 
of anapole or electric dipole moments (fi"7|). 

5 WEAK COUPLING CONSTRAINTS 

From Table |3] one has the matrix element of Ai\ as an expansion in terms of the 
DDH weak couplings. By Eq. (12), one has a corresponding expression for K anapo i e , 
which can be compared to the experimental constraints (Eq. (28)). 

The constraints from 133 Cs and Tl supplement those available from experiments 
in p + p scattering and in nuclear systems. While a larger set of results exists, only 
some of these are generally regarded as being reliably interpretable (|l~4"|). Rather 
precise measurements of the longitudinal analyzing power A z for p + p have been 
made at 13.6 and 45 MeV, and a preliminary result at 221 MeV (where only p 
exchange contributes) is now available. A z has also been measured for p + a at 46 
MeV. There are also two important constraints from nuclei in which observables 
associated with nearly degenerate parity doublets have been measured. In each 
case the nuclear matrix element involved in the mixing has been determined from 
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axial-charge (3 decay (14,46), so that little nuclear structure uncertainty remains. 
The observables are the circular polarization P 7 of the 7-ray emitted in the decay 
of the 1081 keV state in 18 F and angular asymmetry for the decay of the 110 
keV state in polarized 19 F. 
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Figure 8: Constraints on the PNC meson couplings (xlO 7 ) that follow from the 
results in Table |j. The error bands are one standard deviation. The illustrated 
region contains all of the DDH "reasonable ranges" for the indicated parameters. 



A summary of PNC constraints is presented in Table |] and Fig. 8. Although 
the PNC parameter space is six-dimensional, two coupling constant combinations, 
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Table 4: PNC observables and corresponding theoretical predictions, decomposed 
into the designated weak-coupling combinations, with f w = f n — 0.12h p — 0.18/t* 

and h° = h° p + 0.7/i° 



Observable 


Exp.(xl0 7 ) 




h° 


h\ 


h 2 


h° 




AgP(13.6) 


-0.93 ± 0.21 




0.043 


0.043 


0.017 


0.009 


0.039 


Aw (45) 


-1.57 ± 0.23 




0.079 


0.079 


0.032 


0.018 


0.073 


Af (221) 


prelim. 




-0.030 


-0.030 


-0.012 


0.021 




Af*(46) 


-3.34 ± 0.93 


-0.340 


0.140 


0.006 




-0.039 


-0.002 


P 7 ("F) 


1200 ± 3860 


4385 
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-44 


A 7 ( 19 F) 


-740 ± 190 


-94.2 


34.1 


-1.1 




-4.5 


-0.1 


(HAxID/e, Cs 


800 ± 140 


60.7 


-15.8 


3.4 


0.4 


1.0 


6.1 


(\\Ai\\)/e, Tl 


370 ± 390 


-18.0 


3.8 


-1.8 


-0.3 


0.1 


-2.0 



fir-0.12hp-0.18hi and hP-\-0.1h^, dominate the observables, as Table ||| illustrates. 
The la error bands of Fig. 8 are generated from the experimental uncertainties, 
broadened somewhat by allowing uncorrelated variations in the parameters in the 
last four columns of Table || over the DDH broad "reasonable ranges." Note that 
only a fraction of the region allowed by the Seattle Tl constraint is shown: the 
total width of the Tl band is an order of magnitude broader than the width of 
the Cs allowed band, with most of the Tl allowed region lying outside the DDH 
"reasonable ranges" (i.e., in the region of negative f v and positive h° p + 0.7/1°). The 
corresponding Oxford Tl band, which is not illustrated, includes almost all of the 
parameter space in Fig. 8, as well as a substantial region (to the lower left of Fig. 
8) outside the bounds of the figure. 

The weak coupling ranges covered by Fig. 8 correspond roughly to the DDH 
broad "reasonable ranges." Thus the anapole constraints are not inconsistent with 
the most general theory constraints. However in detail, the pattern is disconcert- 
ing. Before the anapole results are included, the indicated solution is a small f n 
and an isoscalar coupling somewhat larger, but consistent with, the DDH best 
value, -(h° p + 0.7h°J^ H ~ 12.7. The anapole results agree poorly with the indi- 
cated solution, as well as with each other. Although the Seattle Tl measurement 
is consistent with zero, it favors a positive anapole moment, while the theory pre- 
diction is decidedly negative, given existing PNC constraints. The Cs result tests 
a combination of PNC couplings quite similar to those measured in A 7 ( 19 F) and 
in AP a , but favors larger values. 

The results plotted in Fig. 8 for 133 Cs are consistent with those of Flambaum 
and Murray (|28|), who extract from the anapole moment an f w about twice the 
DDH best value, f^j?^ ~ 4.6, and point out that theory can accommodate this. 
(The DDH reasonable range is 0-11.4, in units of 10~ 7 .) However, this ignores 
P 7 ( 18 F), a measurement that has been performed by five groups. The resulting 
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constraint is almost devoid of theoretical uncertainty 

-0.6</, -0.11^-0.19/ii < 1.2. (38) 

Allowing hp and to vary throughout their DDH reasonable ranges, one finds 
— 1.0 < /tt < 1.1, clearly ruling out f n ~ 9. Fig. 8 illustrates this, as well as the 
additional tension between Cs, p + a, and A-y( 19 F). 



6 OUTLOOK 



The first conclusion one would draw from Fig. 8 is that additional experimental 
constraints would be helpful. In the case of nuclear experiments, the situation has 
been essentially static for the past 15 years, apart from the p + p A z measurement 
at 221 MeV. However, in the next few years results are expected from experiments 
on the PNC spin rotation of polarized slow neutrons in liquid helium (|49| ) and on 
Ary iun+p — > d+7 (|50"1). Thus we could soon have the information to test whether 
the small- f n solution favored by the p + p and nuclear experiments is correct: this 
is the one consistent solution satisfying the nuclear constraints. 

The anapole situation is less satisfactory: even apart from the question of con- 
sistency with the nuclear constraints, there is some tension between the 133 Cs and 
Seattle Tl allowed regions, though Fig. 8 (which omits the bulk of the Tl allowed 
region) tends to exagerate this disagreement: if the Tl band were enlarged to 2<r, 
it would encompass all of the illustrated region, including the Cs band. While 
the Seattle Tl result favors a sign opposite theory, its broad error bar allows ei- 
ther sign. Taken together, the Seattle and Oxford results are not incompatible 
with any choice of coupling constants with the "reasonable ranges." An improved 
measurement in Tl would clearly be helpful, as would any new anapole measure- 
ment involving an odd-neutron nucleus, which would produce a band in Fig. 8 
roughly perpendicular to that for Cs. There is an atomic PNC effort underway 
using Dy, which has two abundant odd-neutron isotopes, 161 Dy and 163 Dy (30). 
Another possibility may come from a new atomic PNC technique using a single 
trapped Ba + ion: the much larger coherence times and field intensities possible 
with this method compensates for the sensitivity loss stemming from the use of a 
single atom. (The statistical accuracy of traditional methods goes as y/~N, where 
N is the number of atoms.) Ba has two stable odd-neutron isotopes, 135 Ba and 
137 Ba. As with Cs these attempts will have to reach the sub-1% level of precision 
in order to permit an extraction of the differential effects due to nuclear spin de- 
pendence. Alternatively other methods, such as proposals to measure the anapole 
moment directly through E\jM\ interference in hyperfine transitions (51), could 
be developed. 

The underlying issues are not limited to atomic PNC measurements: our under- 
standing of V(e)-A(iV) interactions affects the interpretation of electron-nucleus 
PNC scattering experiments like SAMPLE (|52[), where a similar discrepancy be- 
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tween theory and experiment exists, one that could be associated with an incom- 
plete treatment of PNC effects in the nuclear target. 

Figure 8 shows that the jump in precision achieved in the 133 Cs experiment has 
now provided a constraint on hadronic PNC that is comparable in accuracy to the 
best of the nuclear constraints. This is the reason the lack of overlap in the result- 
ing bands is a concern. The nuclear physics required to analyze the 133 Cs result is 
nontrivial (|5^), and could yet prove to be the source of the discrepancies apparent 
in Fig. 8. The SM calculations described in Section 4 depend on a relation between 
the closure energy and the giant dipole energy that has been tested only in light 
nuclei. It is possible that the systematics for neutron-rich nuclei might be different. 
It is also generally appreciated that spin-dependent operators tend to be quenched 
as model spaces are enlarged to encompass more realistic correlations. Indeed, 
phenomenological single-particle treatments of anapole moments have previously 
envoked phenomenological quenching factors (37). Our SM results are quenched 
relative to single-particle estimates, and would likely be further weakened if those 
calculations could be enlarged. Likewise core polarization results from the RPA 
study of Dmitriev and Telitsin find substantial quenching, relative to single-particle 
estimates (33). All of these results argue that more realistic treatments of correla- 
tions will tend to reduce matrix element values, thus requiring still larger values of 
the weak coupling constants to fit the experimental results. Despite this possibil- 
ity, it is clear that the degree of effort required to complete the 133 Cs experiment 
obligates theorists to invest similar effort in improving their calculations. It is an 
unfortunate situation that the results in Table [|| come from model calculations: it 
is thus difficult to come up with any strategy for quantifying possible errors other 
than the comparisons that can be made when still more ambitious calculations are 
eventually performed. 

This work was support in part by the US Department of Energy and by the 
National Science Foundation. WH thanks the Miller and Guggenheim Foundations 
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